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The choice of a suitable random matrix model of a complex system is very sensitive to the 
nature of its complexity. The statistical spectral analysis of various complex systems requires, 
therefore, a thorough probing of a wide range of random matrix ensembles which is not an easy 
task. It is highly desirable, if possible, to identify a common mathematical structure among all 
the ensembles and analyze it to gain information about the ensemble- properties. Our successful 
search in this direction leads to Calogero Hamiltonian, a one-dimensional quantum Hamiltonian 
with inverse- square interaction, as the common base. This is because both, the eigenvalues of the 
ensembles, and, a general state of Calogero Hamiltonian, evolve in an analogous way for arbitrary 
initial conditions. The varying nature of the complexity is reflected in different form of the evolution 
parameter in each case. A complete investigation of Calogero Hamiltonian can then help us in the 
spectral analysis of complex systems. 

PACS numbers: 05.45+b, 03.65 sq, 05.40+j 

Recent statistical studies in various branches of theoretical physics, ranging from Calogero model of 1-d fermionic 
system [jjj, random matrix (RM) model of disordered systems, matrix models of random surfaces to non- linear sigma 
model of quantum chaotic systems have revealed the presence of a common mathematical structure [^-Q. The 
connecting-web of these various models with each other is well-described in 0. However, so far, the connection of RM 
model with other models was established only for standard Gaussian ensembles (SGE), that is, Gaussian ensembles 
invariant under unitary transformation. This was achieved by showing that distribution of the eigenvalues of the 
ensemble is governed by a Fokker-Planck (F-P) equation HQ similar to that of Dyson's "Brownian" motion model 
H . Through the reduction of F-P equation to the Schrodinger equation, the latter model is already known to be 
connected to Calogero Hamiltonian and thereby to various other models In this paper, we explore RM models 

with non-invariant distributions, and, following the same route as in the case of SGE, connect them to Calogero 
Hamiltonian. This gives us a new technique to analyze the spectral behavior of the quantum operators of complex 
systems. 

The connection between Complex systems and Calogero Hamiltonian seems to be wide-ranging. The eigenvalue 
dynamics of Hermitian operators, for example, Hamiltonians of complex quantum systems e.g. chaotic systems, disor- 
dered systems seems to have an intimate connection with the particle-dynamics of Calogero-Moser (CM) Hamiltonian. 
The latter describes the dynamics of particles interacting via pairwise inverse square interaction and confined to move 
along a real line [jlj, 

here fii is the position of the i particle and V(fii) is the confining potential. Similarly the level-dynamics of the 
unitary operators e.g. time-evolution operator is connected to Calogero-Sutherland (CS) Hamiltonian ||: 

H = - Y £ + ^ Y (>^l) - ?W - 1) (2) 
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where particles are confined to move in a circle thus mimicking the similar confinement of eigenvalues due to unitary 
nature of the operator. The morphological transition caused by the interacting steps on a miscut crystal surface can 
also be modeled by the CS Hamiltonian jTo|| . Here the complexity is thermodynamic in nature. It is already well- 
known that the parametric dispersion of the eigenvalues of the quantum system, with non-integrable classical limit, 
is described by a set of equations similar to the equations of motions of particles, in time, of classical Calogero 
Hamiltonian [ pT[ . This analogy extends also to the statistical properties in the two cases. The parametric-evolution of 
the distribution Pftp(ei, .., £at;t) of the eigenvalues &j of a Hamiltonian H = Hq + tH\ (of size N), with perturbation 
Hi taken from a SGE corresponds to the time-evolution of the distribution Pfrp(rx, ..,r^;t) of positions n's of the 
particles and both the static as well as dynamical correlators of the eigenvalues turn out to be similar to those of 
the particles in CM Hamiltonian (with V(fx) oc fi 2 in eq.(l)) Here j3 refers to the generic symmetry-class of 

the complex systems and therefore the transformation properties of associated RM models (known as GOE, GUE 
and GSE for /3 = 1,2 and 4) p2|,p|. In limit r — > oo, the eigenvalues attain an equilibrium distribution, known 
as Wigner-Dyson, which coincides with the probability distribution of TV-particle coordinates P^p(r;t — > oo) of the 
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ground state of the CM Hamiltonian |^|| . Similar analogies can also be made between evolution of the eigenvalues of 
unitary operators U = J7oe lrM , with M taken from SGE, and CS Hamiltonian [|l3| . This is equivalent to say that the 
stationary and non-stationary states of CSM Hamiltonian correspond to the eigenvalue distribution of the systems 
subjected to random perturbations, strong (t — > oo) and weak (finite r) respectively, and thereby to equilibrium 
and non-equilibrium distribution of SGE. In this paper, we indicate towards a novel connection between the CM and 
RM model: a non-stationary state (finite t) of CM Hamiltonian can also be mapped to the eigenvalue distribution of 
a generalized Gaussian ensemble (GGE); the correspondence is established by identifying a parameter Y for GGE, 
equivalent to time t for CM Hamiltonian. This mapping can then be used to obtain the information about various 
spectral properties of GGEs. 

In recent past, RM ensembles have been quite often used to model the physical systems with complicated interactions 
P,p"2|- The logic which could be given in support of the model is that the missing information about the detailed 
nature of the interactions can be mimicked by randomizing the associated generators of motion, that is, by taking their 
matrix-representations as random matrices. However as the specific details of the complexity of an operator should be 
reflected in the associated RM model, the distribution of the matrix elements can be of various types. For example, 
for a Hamiltonian with chaotic classical limit (least predictability of the long-term dynamics) , the distribution can be 
chosen as Gaussian (the least information ensemble), with distribution parameters to be determined by the associated 
quantum dynamics. The corresponding RM model will thus belong to a generalized Gaussian ensemble with matrix 
elements distribution given by P(H) cx e~fa( H }-fe( H ) (with fi and f% arbitrary functions and H as a typical matrix). 
The SGEs, with matrix elements distribution given by P(H) cx e~ TrH are special cases of GGEs and many of their 
properties are already known. The various features of GGEs have, however, remained unknown so far. The purpose 
of this paper is to suggest a technique to fill in this gap in our information. As for SGE, the nature of matrix elements 
in GGE too depends on the exact symmetry conditions of the Hamiltonian and is again indicated by parameter j3, 
with (3 = 1, 2,4 for a generic matrix element to be real, complex or quaternion \12\. Here we discuss, in detail, the 
properties of GG ensemble of complex hermitian matrices (J3 = 2); the GG ensemble of real-symmetric matrices 
{[3 = l)has been discussed elsewhere |TJ]. We also probe briefly the non-Gaussian ensembles which can serve as good 
models for complex systems with various conditions on the associated quantum dynamics. 

We proceed as follows. Our technique is based on the statistical evolution of the eigenvalues of a GG ensemble with 
respect to a change in their distribution parameters. This requires a prior information about the effect of a small 
change in the matrix element on eigenvalues and eigenvectors; the related study is given in section LA. These results 
are then used to obtain, as described in section LB, the distribution of eigenvalues P (/it, V) of a matrix H taken from 
a Gaussian ensemble, non-invariant under unitary transformation. The evolution of the eigenvalues is governed by a 
partial differential equation which, after certain parametric redefinitions, turns out to be formally the same as the F-P 
equation for the Brownian motion of particles in Wigner-Dyson (WD) gas jl2| . The section II contains the details of 
the reduction of the F-P equation to the Schrodinger equation of CM Hamiltonian and a mapping of their respective 
correlators. The section III deals with the application of our technique to some important physical processes e.g. 
localization and a brief discussion of our technique applied to a few other important matrix ensembles is given in 
section IV.. We conclude in section VI which is followed by the appendices containing the proofs of some of the results 
given in the main text of the paper. 



I. EIGENVALUE DISTRIBUTION OF GENERALIZED GAUSSIAN ENSEMBLES 

A. The Change of Eigenvalues and Eigenfunctions 

The eigenvalue equation of a complex Hermitian matrix H is given by HU = UA with A as the matrix of eigenvalues 
A n and U as the eigenvector matrix, unitary in nature. As obvious, a slight variation of the matrix elements of H 
will, in general, lead to variation of both the eigenvalues as well as the eigenvectors and associated rates of change 
can be obtained as follows; 

As A„ = U n iHijU*p the rate of change of A n with respect to Hki- S (with s referring to real, s = 1, and 
imaginary, s = 2, parts of Hu) can be given 

f)\ ? s_1 

[u ln u* kn -(-iyu? n u kn ]. (3) 



dHki-s gn 

where gki = 1 + Ski ■ This can further be used to obtain the following relations (Appendix A) 



E E JWT-" 1 ' 1 - = E HuUmUL = A„ (4) 

k<l s=l ' k,l 
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and 

k<l s=l 1 

For our analysis later, we also require the information about the second order change of an eigenvalue with respect 
to a matrix element and, therefore, the rate of change of one of the eigenvector components with respect to Hki- This 
is given as follows (Appendix B), 

f)TJ AS — 1 1 

7uf^ = it E a — ^-u pm {u* km u ln + (-iy +1 u; m u kn ) (6) 

On ki- s g k i *—f A„ — A m 
and now by using eqs.(3,6), One can show that (Appendix C) 

E^E^^Er^x- (7) 

k<l s—1 rn 

For the real-symmetric case, the corresponding relations can be obtained by using U + — U T (as eigenvector matrix 
is now orthogonal) in eqs.(3-7) and taking Hij-2 = for all values of i,j (see f|). 



B. The Evolution Equation For the Eigenvalues 

Let us consider an ensemble of complex Hermitian matrices H, with matrix elements Hki = Hki ; i + iHki-^i^- — $ki) 
distributed as Gaussians with arbitrary variances and mean-values; the variances of real and imaginary parts of a 
single matrix element also need not be same. Thus we choose the distribution p{H) of matrix H to be following: 

2 

p(H, y, b) = Ccxp(- J2 E a M;s{H U ; S - b kl , s f) (8) 

s=l k<l 



with C = Ilfc<i y t?" as ^ e normalization constant, y as the set of the coefficients y k i;s = Oikl;s9kl 

and b as the set of all bki- s - Note that such a choice leads to a non-random complex Hamiltonian (Hki = bki-i + ibki-,2) 
in limit aki-i, cxm-2 —> oo and therefore can model various real physical situations such as switching of disorder in a 
non-random Hamiltonian e.g. metal-insulator transitions. 

Let P(p, y, b) be the probability of finding eigenvalues Xi of H between fii and p,{ + dpi at a given y and b, 

N 

P(H, y,b)= / J] S(jH ~ h)p{H, y, b)dH (9) 

i— 1 

As the a-dependence of P in eq.(9) enters only through p(H) and a ®£ = [(2afc; ;s ) _1 — (Hki- S — bki- s ) 2 ] p = 

(2au-,s)~ P + {Hkl-s — bkl-sJ-ffH 

as follows 



dp 



with g ^ p = — d ® p , a derivative of P with respect to aki- s can be written 



'"' - • — f l\S(^ - X^Hki-.s^^dH + fT[8(j ii -X i )b k i ;B -^-dH 

ki-s J XJ r oHki-s 2aki;s J AJ r ob k i- s 

l—l 1=1 



dahlia 20tkl:a %&kl 



(10) 



The second integral in eq.(10) is equal to bki_ s gg^- ■ The first integral can also be simplified by using integration 
by parts followed by a use of the equality g Hkl , s = - L„=i a M „ 31^77 : 

/ TT 5(jH - X i )g M ..-^-dH = - / dlUi %i - Ai) gfc dg _ fTJ 5 ( Jii _ X .) pdH = J p (11) 

i=l ' 1 
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where 



N fi r N f>\ 

n=l i=l 



Substitution of eq.(ll) in eq.(10) then gives 



dP dP 



' da k i-s 



'db 



kl;s 



(12) 



(13) 



Our aim is to find a function Y of the coefficients a/y ;s 's and i>jy ;s 's such that the evolution of P(u, Y) in terms 
of Y satisfies a F-P equation similar to that of Dyson's Brownian motion model (Wigner-Dyson gas) p,^2[|. For this 
purpose, we consider the sum 2 X)fc<2 (7 ~ 9kl&kl;s) &kl;8 where 7 is an arbitrary parameter and thereby obtain 

following relation 



yu-s) 



s=l k<l 



dP dP 

2ykl;s— Okl;. 



E 



k<l 



k<l 



' dyu-s ' db k i. s 

As shown in Appendix D, the first term on the right hand side of eq.(14) can further be simplified, 

s=l k<l n 

The second term can similarly be rewritten as follows (Appendix E): 



(14) 



(15) 



E vki;shi, s - - E 

s k<l n " n 



— v 



f_l 7n f_l 7] 



P - },ykl;shl;6 



k<l 



dP 
db k i. s 



(16) 



where (3 — 2. Using both the equalities (15) and (16) in eq.(14), we obtain the desired F-P equation 



dP d d 

w = 7E^^» p )+E^ 







\ ~" 



p 



(17) 



Here the left hand side of above equation, summing over all yki-, B and bu-s, has been rewritten as §jy with Y given 
by the condition that 



dP , .BP BP 

s k<l y ,s s k<l 



dY 



' dbki.s 



(18) 



By using the orthogonality of eigenvectors and following the same steps, it can be proved for real-symmetric case too 
(now (3—1) It is worth noting that the eq.(17) is same as the evolution equation for the eigenvalues of Brownian 
ensembles. It is also similar to the one governing the transitions between any two universality classes of SGE caused 
by a random perturbation of strength r (with t — > Y) fipl . 



C. How to Obtain the Complexity parameter Y: 



The variable Y, a function of relative values of the coefficients ct^s's and 6fe;- s 's, is a measure of the degree and 
nature of the complexity of a system and can therefore be referred as the "complexity parameter". For the case 
discussed here (eq.(18)), Y can be obtained by the following method. 

We define M = 2N 2 variables (Yi, ..Ym) as the functions of all t//y ;s 's and bki-, s 's such that the condition given by 
eq.(18) (where Y = Y\) is satisfied. This is indeed possible by using the orthogonal (Jacobi) coordinate transformation 
between variables {5^}i=i,..,M and {j/w ;s , bki-, s }k<l;k,l=i,..,N;s=i,2 defined by following rule, 
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M 



for 



M 



(19) 



where X« 



±ln T 



of integration. 



Here coefficients a. 



Cj for j < N 2 and Xj 



must satisfy the relation ^ 



Cj for j > N with e, as arbitrary constants 

*y muo. a »u«,ij uuv, i^.uuiuii z-fj=i a u = <^i which is a necessary condition for 
orthogonality but not sufficient to get the right form for ^7. With D being the functional derivative of Y^s with 
respect to Xj's, we also need the elements Dy- of its inverse to be unity. One way to achieve this is to set all adjuncts 
of the matrix elements Jy- equal. Now by choosing ay also equal, ay = M~ l , we are left with M conditions for , 
i ^ 1, which can easily be fulfilled. 

The form of Y — J2j a ijXj, fulfilling condition (18), can therefore be given as 



Y = 



2N 2 



EE 

k<l s=l 



iln 

2 \Vki 



Vku. 



■1\ 



~\n\b k i;s\ 
7 



C 



(20) 



with C = M _1 £\ c^. 

As obvious, this method is applicable only for the case when the prefactor associated with a derivative of P with 
respect to a variable r in eq.(18) depends only on r (r can be any one of the yki- s or bki- s )- Our studies on the ensembles 
more complicated than eq.(8) show that the prefactors can also depend on a combination of various r variables. This 
requires a more general method to obtain Y which can also be used for the case discussed here (Appendix F). 



D. Determination of P(/j,,Y) 



The eq.(17) describes an evolution of the eigenvalues of GGE due to changing distribution parameters of the 
ensemble which can be solved, in principle, to obtain P(/i, Y) for arbitrarily chosen initial values of the parameters. 
If the ensemble corresponding to initial set of the parameters is referred as Hq, an integration over Hq would lead to 
P(fJ,, Y), free of initial conditions. In fact, it can be shown that 



P(p;Y) = (4nY)- N2 /* 



exp 



4F 



tr(fx - U+iMiUf 



/MIAMI 00 d^dU 



(21) 



where [1q is the set of eigenvalues of the initial matrix Hq, with (3q given by its symmetry conditions, and U is the 
integral over unitary (or orthogonal) space of matrices. 

To show that eq.(21) is indeed a solution of eq.(17), we study a general case. Consider a partial differential equation 
for a function F(A; t) defined in the matrix space of N x N Hermitian matrices A 



^=[ViF + V.(AF)] 



where 



V 2 4 



d 2 



IE 



2 ^ dA% 



and 



v-(^)=E^(^) 



i<3 



(22) 
(23) 



with the initial condition F(A; 0) = f(A). This equation is known to have a unique solution (see page 174 of |jT2f ), 

F(A;t) = J K(A,B,t)f(B)dB (24) 

(25) 



where 



K(A,B,t) = (4wt)- N /2 exp 



tr(A-B) 2 

At y ' 



where B is a N x N hermitian matrix. Depending on the nature of both A and B, we can choose a special class 
of eigenvector matrices U a and Ub (for A and B real-symmetric, complex hermitian or symplectic, Ua and Ub are 
orthogonal, unitary and symplectic matrices respectively) such that 



A = U s A aU A and B = U s B bU B 



(26) 



where a = [a»<Sjj], b = [biSij] are diagonal matrices with ai and bi as the eigenvalues of A and B respectively and 
U s = U + or U T or U R depending on whether U is an eigenvector matrix for a complex Hermitian, real symmetric 
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or symplectic matrix ]l2]| . Let 0a and 0b give the number of components of a typical matrix elements in A and 
B respectively. Changing the variables from matrix elements to the N eigenvalues and f3N(N — l)/2 angle (i.e 
eigenvector) parameters on which Ub depends, we have 



AB = |A(6)| /3t db dU B with db = ]Jdb z and A(b) = Y[(b z - b 3 

i<3 



1=1 



The substitution of these relations in eq.(24) gives us 

1 



F(A; t) = (2nt)- N I 2 



exp 



2t 



tr(a-U s bUy 



/(&,C/ B )|A(&)| ft d&df/i 



(27) 



(28) 



where U = UbU s a and II s = UaU b . Now if /(&, Ub) is independent of Ub then F(A; t) is also independent of Ua- 
This helps us to rewrite the eq.(28) as follows, 



F{a;t) = {^ty N2 ' 2 / G(a,b,t)f(b)\A(b)f»db 



(29) 



where 



G(a, b,t) — J exp 



4f 



tr(a - U s bUf 



dU 



Here the integral is over the group U of orthogonal, unitary and symplectic matrices respectively. Further the 
Laplacian V\ can also be written in terms of eigenvalues and angle parameters of A (see appendix A. 5 of | fl2"|| ) 



V^) = ^£^|A<»)|-£ + V 



2 

Ua- 



(30) 



By the substitution of eq.(30) in eq.(22) and using independence of F(a; t) of Ua, one can rewrite eq.(22) as follows, 



dF{a-t) 



|A(o)|A 



E 



da,; 



(31) 



with P(a; i) given by eq.(29). Now by using the equality J2i ^l^( a )l /3a = 0> eq.(31) can be reduced in the following 
form: 



OF ^ d . _. ^ a 



9a, 



9a ; ^-f a,- - cii 



F 



(32) 



which is similar to eq.(17) with — > /i;, t — > Y", 7 = 1 and F —> P. The joint probability density P can therefore 
be obtained by evaluating the integral (29). However, so far, the integration could be performed only for the unitary 
group of matrices • 



E. Steady State, Level Density and Correlations 

The steady state of eq.(17), P{p, 00) = P^ — \A(p)\ /3 e~^ corresponds to Y — Y a — > 00 (with Y a as the 

complexity parameter of initial ensemble) which can be achieved by two ways (for finite Yq values). The first is when 
almost all yua ~ > 7 and yui-,2 — * 00 (for finite bki-i and bki-2 values ) which results in a GOE steady state. The 
second is when almost all yui-.i — > 7, Vkl-,2 - > 7, resulting in a GUE. This indicates that, in the steady state limit, 
system tends to belong to one of the SGEs. The eq.(17) can, therefore, describe a transition from a given initial 
ensemble (with Y = Yq) to either GOE or GUE with Y — Yq as the transition parameter. The non-equilibrium 
states of this transition, given by non-zero finite values of Y — Yq, are various Gaussian ensembles corresponding to 
varying values of the coefficients yki ;s and bki- s - For example, the choice of the initial ensemble as GOE (almost all 
yu-,i = liVM;2 ~ ^ 00 initially) and a decrease of yki-,2 (from oo — > 7 while keeping yki-i fixed) leads to GOE — > GUE 
transition with intermediate ensembles as those of complex Hcrmitian matrices. Similarly Poisson — * GUE transition 



G 



can be brought about by choice of the initial ensemble as Poisson (almost all yki;i,Vkl;2 —> oo for k ^ I, yuk;i = 7j 
Vkk-2 = 7 and bki- s — for all k, I, s values) and by varying both yui-i and yu-,2 upto 7. As clear from above, 7 fixes 
the variance of the final ensemble and an arbitrariness in 7 leaves the latter arbitrary. This however does not affect 
the statistical properties of the intermediate ensembles. 

The eq.(21) for P{fi,Y) can be used to obtain n th order density correlator ../x n ; Y), defined by R n = 

(jv-n)! / P{^^ ^)d/in+i"d/ijv- {Rn can also be expressed in the form < v{fJ,i, Y)..v{(i n , Y) > with v{(j,, Y) = 
A -1 J7(5(/i — jtij) as the density of eigenvalues and < .. > implying the ensemble average). Here note that the 
analogy of eq.(17) to that of Dyson's Brownian ensembles implies same form of P for both the cases and therefore R n . 
A lot of information already being available about level-density and various correlation for Brownian ensembles, it can 
directly be used for ensemble described by eq.(8). Thus, as for BE, a direct integration of F-P equation (17) leads to the 
BBGKY hierarchic relations among the unfolded local correlators Rn{ri, ..,r„; A) = LimA — > 00 ^^^y)']] Ri(p^- y ) 
with r = f Rx(fj,;Y)dfJ. and A = (Y - Y )/D 2 {D = R^ 1 ; the mean level spacing) fy|, 

9Rn _ \ - d 2 R n _ aST^ ( \ _ aST _iL f°° Rn+l /oo% 

OA ^ dr 2 P 2-, Or, { Tj rj '2- dr, J x r, - r, 

(here for simplification, 7 is chosen to be unity). As can be seen from the above equation, the transition for R n 
occurs on the scales determined by Y « D 2 and a smooth transition can be brought only in terms of the parameter 
A, obtained by rescaling Y by D 2 . On the other hand, for the corresponding scale is given by Y w N D 2 . This 
implies, therefore, during the transition in i?„, the density R\ remains nearly unchanged; this fact is very helpful in 
unfolding the correlators R n . For n = 1 and in large A- limit, above equation reduces to the Dyson-Pastur equation 
[p| for the level density < v{pi,Y) >= N~ 1 Ri 

o<»M>_ a d / N -„ [^<vW)>\ „ !/;r ( , 4) 




which results in a semi-circular form for v\ v{r) = -^-{v 2 — r 2 ) 1 / 2 with rj 2 — 4A(1 + Y 2 ) The application of 



super-symmetry (SUSY) technique fllq] to ensemble (8) gave a similar result (also see section 4.3 of 12 ) 



II. CONNECTION TO CALOGERO HAMILTONIAN 

A similarity transformation followed by a Wick rotation converts the F-P equation into a self-adjoint form ||. This 
can be seen as follows. The F-P equation, in general, can be expressed in a form 

where P is a F-P operator with non-negative eigenvalues. Here \Py > is a general state of operator P at "time" 
Y and its projection in eigenvalue space can be obtained by the usual operation P{fj,,Y) =< fi\Py > (with \x as 
set of the eigenvalues). Let P{fi, Yo) =< fJ-\Po > be the equilibrium probability. One can further define a vector 
< 0| = J d/i < satisfying < 0\P = thus implying the conservation of probability in "time" Y in this state (the 
ground state). The F-P operator can now be hermiticized through a similarity transformation S~ 1 PS = H where S is 
Hermitian and invertible operator depending only on the eigenvalues. Thus the ground state condition must be given 
by HS\0 >= (as P + \Q >= 0). Let the effect of similarity transformation on the state \Py > and \Pq > is expressed 
by \ip >= S' _1 |Py > and \ipo >— S' _1 |Po > respectively. The similarity transformation of eq.(35) will then give the 
desired form d 'J^ = —H\tp >; the ground state \ipo > must now satisfy the condition H\ip >= 0. The comparison 
of the two different forms of the ground state condition gives \ipo >= S\0 > and therefore \Pq >= S 2 \0 >. 

In the case of F-P equation (17), H turns out to be CM Hamiltonian (eq.(l) with r"j — > /U,) and has well-defined 
eigenstates and eigenvalues 0Jl^l. As well-known, the particles in CM system undergo an integrable dynamics, thus 
implying a similar motion for the eigenvalues too. Here H being a generic member of GGE, this result is valid for all 
systems with interactions complicated enough to be modeled by GGE. 

The "state" ip or P(/i, Y\Hq) can be expressed as a sum over the eigenvalues and eigenf unctions which on integration 
over the initial ensemble Hq leads to the joint probability distribution P(/i, Y) and thereby static (at a single parameter 
value) density correlations R n . The above correspondence can also be used to map the multi-parametric correlations 
to multi-time correlations of the of CM Hamiltonian. For example, the parametric correlation < Q a {Y)Qb{0) >, 
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for a classical variable Q(Y) with [Q, S] = can be mapped to the corresponding ground state correlation of CM 
hamiltonain < ifio\Qa(Y)e~ YH Qb(0)|V>o >■ This follows because 

< Qa(Y)Q b (0) >= J QaQbPQi, F)d/i = J < ti\Q a Q b \P Y > d M (36) 

now as the evolution of \Py > with respect to Y is given by \Py >= Se H S^IPq >, one has 

< Qa(Y)Q b (0) >=< 0\Q a Se' YH S' 1 Q b \P o >=< MQa^ YH Q b \^ > (37) 



III. APPLICATION TO PHYSICAL PROBLEMS 



The given ensemble (8), referred here as "G", is represented by a point Y in the parametric-space consisting of 
distribution parameters and various transition curves may pass through this point. The question therefore arises 
which curve should be chosen for the studies of the properties of G? The answer is the one which does not leave 
any arbitrariness behind and if there are more than one such curve, each one of them should give same answer for 
various fluctuation measures of G. This criteria for the right choice are based on the symmetry properties of ensemble 
G, that is, the nature of all au and b^i with end-points (the final and initial ensemble, referred here as "F" and 
" O" respectively) chosen in such a way that the values corresponding to G occur during the variation of distribution 
parameters from one end to the other. Further the chosen transition should preferably be the one whose properties 
are already known and can therefore tell us about G. For many GGE described by eq.(8), above criteria is satisfied by 
choosing F as a SGE with variance < >= 2 < F 2 - > = 7~\ 7 < mm{y k i- s [G)} , k, I = 1, 2, .., AT, s — 1, 2, and O as 
an ensemble with each a k i [O] given by the maximum value taken by the functional form of the corresponding a k i [G] . 
However, as explained in following examples, O can also be chosen as some other ensemble. For example, If G is an 
ensemble of real-symmetric matrices H represented by p{H) cx exp[— X)/c<; a kiH k i\ with finite but different values for 
all aw, the Poisson — > GOE curve is more suitable for its study rather than GOE — > GUE. Here the GOE ensemble 
is described by < >= 2 < Ff, >= -f^ 1 with 7 as the minimum value among given y k i[G]s. However if various a k i 
in the above example can also take complex values, the ensemble can now be chosen on any one of the curves, namely, 
Poisson -> GUE or GOE -> GUE. Here now GUE can be chosen as < F* >= 2 < F 2 ;1 >= 2 < F? y2 > = 7- 1 . 
The GOE for the second curve can be chosen as the one with < >= 2 < 0?-^ >= q^ 1 and < 0^. 2 >= with 
q = max{yjj;i [G]}. Similarly, for Poisson — > GUE curve, the initial ensemble may be taken as one with < 0\ >= 7 _1 
(or q _1 ) and < 0\^. x >=< Ofj. 2 >= for i 7^ j. The reason for the choice of the two transitions is due to availability 
of the results for their 2-point correlation R 2 [|l3| : 

For Poisson — > GUE 



R 2 (r;A) - i? 2 (r;oo) = -J dx J dz cos(2vrrx) exp [-8tt 2 Aj:(1 + x + 2zy/x)] (^~^ 



v/(l-z 2 )(l + 2zjx) 



2z^/x 



(38) 



and for GOE -> GUE 



R 2 (r; A) - R 2 (r; 00) = -\ f dx [ dz x sm(rx) exp [2A(x 2 - y 2 )} (39) 

ti" Jo h V 

where R 2 (r, 00) = 1 - si ^ r) (the GUE limit). 

It is obvious therefore that if Ai and A 2 are the parameter values for the ensemble "G" on Poisson — > GUE and 
GOE — > GUE curves respectively, one should have R2,p^u(r; A = Ai) = i?2,o— >!/('"; A = A2). This would require 
an intersection of two curves in the R 2 — A space which however is possible. This is because the GOE can occur 
as an intermediate point in Poisson — > GUE transition. The GOE — * GUE curve can also appear as a part of the 
Poisson — > GUE curve; thus the choice of two different initial ensembles here corresponds only to two different origins 
of dynamics on the same curve. 

The parameter 7, which determines Y as well as the variances of F 1 , enters in calculation at step given by eq.(14) and 
can be chosen arbitrarily. As suggested by eq.(17), the choice of different 7- values corresponds to different y-values 
as well as the transition curves with end-points of same nature but different variances; this, however, would not imply 
different properties for the ensemble G (Appendix G). Similarly the F-P equation is although independent of the 
choice of the initial ensemble, the latter is required for determination of the correlations of G. The possibility of an 
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arbitrary choice of O may seem to imply a certain arbitrariness left in the correlation of G. However the choice of 
two different initial ensembles corresponds only to the two different origins of the dynamics approaching to the same 
point in the parametric space. 

It will be clarified by the examples given below. 



A. Anderson Transition 



Using above method, the transition parameter for a metal-insulator transition as a result of increasing disorder 
can exactly be calculated. To see this, let us consider the case of a d-dimensional disordered lattice, of size L, in 
tight-binding approximation. Here, in the configration space representation of the Hamiltonian, a N x N matrix of 
size N — L d , the diagonal matrix elements will be site-energies e,-. The hopping is generally assumed to connect only 
the z nearest-neighbors with amplitude t so that the electron kinetic energy spread or bandwidth is zt. This therefore 
results in sparse form of the matrix H . We first consider the case of L — > D transition brought about by decreasing 
diagonal disorder only. In this case, site-energies e; are taken to be independent random variables with probability- 
density p{ej). In the Anderson model of metal- insulator (MI) transition, p(e) was taken to be a constant W~ x 
between — W/2 to W/2. Various physical arguments and approximations used in this case led to conclusion that here 
all the states are localized for W > AKt\n{Wr) with K as a function of z and d. 

However, as well-known now, MI transition does not depend on nature of pie) and latter can also be chosen as 
Gaussian; the type of p(e) affects only the critical point of the transition. The p(H), for any intermediate state of MI 
transition brought about by diagonal disorder, can therefore be chosen as in eq.(8) with au — > oo, bki = — t for k I, 

afefc = « and bkk — for all k which results in Y = T m |2a-7| — 7 _1 Adni + C. Here K is total number of 

the sites connected and depends on the dimensionality d of the system. The system can initially be considered in an 
insulator regime where all the eigenvectors become localized on individual sites of the lattice (strong disorder limit). 
This results in a diagonal form of the matrix H with the eigenvalues independent from each other. The insulator 
limit can therefore be modeled by ensemble (8) with au — > oo for k ^ I, akk = otQ (for all k- values) and bki — * (for 
all k, I), giving, Y"o = ijv m | 2a°- 7 | + ^ ( as K=0 in the insulator regime). The decrease of the diagonal disorder, that 
is, an increase of akk from an to some finite values (while ctki, k ^ I, remains infinite throughout the transition) will 
ultimately lead to metal regime with fully delocalized wavefunctions. The eigenvalue distribution of H in the regime 
can be well-modeled by the SGE; let it be described by um(> an). Thus for the study of transition in this case we 
should choose 7 = 2a m- The transition parameter can now be given as follows, with the mean level spacing D oc -^=, 



D 2 4 



a\a - a M \ K 

m— \ r — ^ m * 

ao\a-aM\ Nam 



(40) 



As obvious from the above, the transition is governed by relative values of the disorder and the hopping. Here A — ► 
leads to fully localized regime which corresponds to following condition on a and t 

, a a — an K , . 

In— + °- = - hit (41) 

a a M Na M 

The eq.(40) gives, therefore, the condition for the critical region or mobility edge (|j> — > finite as N — > 00). As 
\ a ~ a °\ << 1 even for large a- values, the condition is always satisfied if — > 0. This explains the localization of 
all the states in infinitely long wires (or strictly 1-d systems where K « N) even for very weak disorder. With 
increasing dimensionality d, connectivity K of the lattice and thereby the possibility of |A| >> and the delocalized 
states increases. The A can similarly be calculated when off-diagonal disorder is also present. 



B. 1-D, Quasi 1-D, Periodic 1-D disordered and Chaotic systems 



In 1-D geometry of a solid state system e.g a chain of N interacting sites, in tight binding approximation, the 
long-range random hopping leads to a banded structure of the matrix, known as random banded matrix (RBM) 
pl| , p2| . Here the effectively non-zero, randomly distributed, matrix elements are confined within a band with its 
width governed by the range of hopping. The 1-D periodic geometry e.g. a chain of interacting sites joined into 
a ring leads to a periodic RBM in which all non-zero matrix elements belong to three regions: a band along the 
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main diagonal, the upper right corner and the lower left one |pl[ . A real disordered wire has finite cross-section 
(referred as quasi 1-D geometry) and therefore allows for propagating modes with different transverse quantization 
numbers frequently referred as transverse channels. This case can again be modeled by RBMs with band-width given 
by number of transverse channels 0,E2]. In the case of dynamical systems too, exhibiting strong chaos in classical 
limit, the generic structure of Hamiltonian matrix in some basis is banded and matrix elements can be assumed to be 
pseudo-random pi|| . For example, the Hamiltonian of quantum kicked rotor turns out to be a RBM in momentum 
basis |||. 

In all these cases, nature of the disorder or associated randomness decides the nature of the distribution of the 
matrix elements. The physical properties of such systems can therefore be analyzed by studying the distribution of 
the eigenvalues of associated RBMs. The most studied type of RBM is that with the zero mean value of all matrix 
elements and variance given by < i?„ m >= v 2 a(\n — m\/b) where a(r) is some function satisfying the condition 
lim r ^ 00 a(r) = and determines the shape of the band pSq] . For large but finite size of the matrix N >> b >> 1, 
its statistical properties were shown (by SUSY method) to be determined by the scaling parameter b 2 /N with the 
transition parameter scaling as Nf(j^) pl| . 

The transition parameter for the RBM can also be calculated by our method. Let us first consider the simplest 
case i.e. Rosenzweig-Porter Model where all the off-diagonal matrix elements are distributed with same variance 
which is different from the diagonal ones. Let us take onj;i=£j[G] — 2(1 + fx) and a^[G] = 1 with fi > 0; thus 
mm{yij [G] } — 2 and we can choose 7 = 2. This GGE can therefore be mapped to a Brownian ensemble, with 
Y — Y"o = — ^^y 1 ln| 1 — yq^jl ~ jiz for \i > 1, appearing in a Poisson — * GOE transition where the initial matrix 

EH 2 
. . i " and the final, stationary state, obtained for large A-values, is 

P(H) cx e~^ TlH2 . Now as i?i « §|J|, the D 2 « l/N and therefore A w ^ which implies that the GGE will 
have an eigenvalue statistics very different from that of Poisson or GOE only if \i « cN (c a finite constant). For 
fi > cN, A — > and for /1 < cN, A — » 00 for N — > 00 and thus the GGE behaves like a Poisson ensemble in the 
first case and like a GOE in the second; this result is in agreement with the one obtained, in pjj, by using NLSM 
technique. (Note in ref. p5|| , D is taken as D cx l/N, which gives A«|- and therefore GOE and Poisson ensemble 
result for fi < cN 2 and [i > cN 2 respectively). 



Consider the ensemble with exponential decay of the variances away from the diagonal i.e au = e' fe l ^ b ,k < 
1,1 << b << N. Thus, again 7 — 2 and the final ensemble is a SGE with P(H) = c -i TTH ~ and therefore Y = 
~aW i ^2i<j=i — ^yij 1 ^ ^ ' Here the initial ensemble is that of the diagonal matrices with a Poisson distribution 
of the eigenvalues which corresponds to ya[0]^= 2 and 2/»j;iyy[0] — > 00 (this being maximum value of yw[G]) giving 
Y» = -jm Eili Ml ~ 7^ X [0]| + C. Thus Y - Y = - r)ln|l - 2e~ r / b \ w A. As R x « y/N, the 

transition parameter for infinite system (N — > 00) turn out to be A = w b (see ]14|| ) which reconfirms that, in 
infinite systems, the transition is governed only by the band- width b ]2^ , |2l| . 

Another case of importance is the ensemble with power law decay of variances Hij — Gijd(\i — j\) with G a 



typical member of SGE (< Gf 



>= 2 < G 2 j 



>= v 2 ) and a(r) 



respectively (known as PRBM model with P stands for power) 



= 1 and (b/r) a for r < b and 
26). This corresponds to yij ■ 



therefore 7 = min{j/y} = 
the same. Now as = y, 

A = D- 2 (Y-Y ) 



> b (b » 1) 
v*a*f\i-j\) and 



Again as for the exponential case, the choice of initial and the final ensemble remains 
7(0 2<T (with r = \i- j\), we get 



1 N ( 

-iiv E(^-> 1 

r=b+l \ 

4^iW hi. 

A 2^1 A 



(42) 



3=1 J 



b/N 
1 



dx(l — x)x 



-2ja 



2(l-2ja)(l-ja) \N 



2ja 



1 



1 



(l-2ja)N 2(1 -ja)\N 



(43) 



Thus, for large A-values and a < 1/2, A(oc N 1 ~ 2cr ) is sufficiently large and the eigenvalue statistics approaches SG 
limit. At a — 1/2, the statistics is governed by the parameter b 2 /N instead of N only. For a = 1, the non-zero, finite 
A-value (A cx b even when N — > 00) leads to an eigenvalue statistics intermediate between that of SGE or Poisson. 
For a > 3/2 with N — > 00, A — > therefore, the eigenvalue statistics approaches Poisson limit, A being very small. 
All these results are in agreement with those obtained in p6l by SUSY technique. 
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Another type of RBMs often encountered in atomic and nuclear systems are those with the non-zero mean value 
of all matrix elements and with variance given by < H 2 m >= v 2 a(\n — m\/b); the transition parameter for them can 
also be obtained as for the above cases p?| , |28| , ^3| , p^ , ^0| . 



C. Quantum Hall Case 



A Quantum Hall system without disorder has all the states degenerate within a single Landau level. The introduction 
of disorder leads to a broadening of the levels (also termed as diagonal disorder) as well as random hopping between 
them (off-diagonal disorder) and a competition between the two leads to a L — » D transition. Note this is different 
from the Anderson model where the L— * D transition is caused by the competition between diagonal disorder and 
non-random hopping (bandwidth) . The Y x Y Hamiltonian matrix in presence of disorder therefore belongs to 
an ensemble far more complicated than eq.(8), known as random Landau matrix, as now various matrix elements are 
no longer independently distributed: p{H,y,b) = Cexp[- J2s=i J2k : i-,k<i H k i- s { a ki;sH k i- s - V^.,.,. , &y« ;s -ffy>)] with 
C as the normalization constant and y and b as the sets of inverse of variances yki- s = &ki;s9kl & n d coefficients bijki-s 
respectively with gu = 1 + Ski- Here Ylij ^ijfcZjs will imply that the summation is over all possible pairs of indices 
{i,j} such that the pair {i,j} ^ {k, 1} or {I, k} [pl[ . In this case too, one can show that the eigenvalue distribution P 
satisfies eq.(17) but the condition for the determination of Y is no longer given by eq.(18); the details will be presented 
elsewhere. 



D. Critical Ensemble and Multifractality of Eigenvectors 

Recent studies of some metal-insulator transitions revealed that the energy level statistics in the critical region 
is universal and different from both Wigner-Dyson as well as Poisson statistics. The eigenfunctions associated with 
the critical statistics show multifractal characteristics [[32-34 . The level number varaince £ 2 (Y) is believed to be 



an important indicator of this critical behaviour with its asymptotic linearity in the mean number of levels TV [p5| 
£ 2 (Y) =< (SN) 2 >— xN, x < 1- The critical statistics, therefore, governs the spectral fluctuations that are weaker 
than for the Poisson statistics (£ 2 (Y) = Y) but much stronger than for the Wigner-Dyson statistics, (£ 2 (Y) = lnY). 
Later on remarkable similarities were found between the spectral statistics of a number of dynamical systems e.g 
pseudointegrable billiards and the critical statistics near the mobility edge |3(| . However such a critical region being 
inaccessible either perturbatively or semiclassically, a different tool was required to probe into it. This led to the 
suggestion of a RM modelling of this region J34|. The Y x Y matrices in this model are Hermitian and matrix 
elements are Gaussian distributed with zero mean and the variance given by 



< (Hij) 2 >-- 



■3 



B 



(44) 



Using SUSY technique, it has been shown pq] that for large B-values {B >> 1), this ensemble behaves like a SGE for 
a < 1 and as a Poisson for a > 1. The case a = 1 is believed to be of special relevance as it supports critical statistics 
and multifractal eigenstates; the application of SUSY technique gives -R2M ~ 1 — j^gs S i n h°(ir/4B) an< ^ ^ 2 (^0 ~ xN 

©mi; 

The existence of the ensembles with critical statistics is indicated by our technique too. The Y-dependence of 
the transition parameter A, entering through Y and the mean level-spacing D, causes the transition to reach the 
equilibrium in limit N — > 00 for finite, non-zero Y-values. In some cases, however, the Y-dependence of Y may 
be such that it balances the one due to D, thus resulting in an Y-independent A (as shown in section III.A,B) 
and therefore critical statistics. As can be seen from eq.(20), A for the ensemble, given by eq.(44), is also Y- 
independent for a = 1; here the ensemble appears as an intermediate point between Poisson — > GUE transition with 

Y — Yo = J2^=i(^ ~~ r ) m (l + (r) 2a ) an d ^ behaves as in the case of PRBM model discussed above, showing 
criticality for a = 1. The correlation R 2 for the ensemble (44) can therefore be given by eq.(38) which for large 
A- values (for all r), can be approximated as follows p0|j39||: 



R 2 {r,A) = 1 + — ^ + -4-^[cos(27rr)e- 2 T - 1] (45) 

7r 2 A 2 + r 2 2ir 2 r 2 L 



= 1 + -5T + 7r4^[e" 2 ^ - 2e- 2 T S in 2 (7rr) - 1] (46) 
7T 2 A 27r 2 r 2 
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1 sin 2 7rr 6 sin 2 (7rr) . r-. 

~ 1 + — tt" ~ 1 tt o , /- , ; for r « V A) (47 

tt 2 A 7r 2 r 2 e 2^ 7T 2 A s inh 2 (rV6/A) ; V ; 

which is similar to the result given by SUSY technique. However, for A >> r » VA~, our method gives 1 — R%(r, A) = 
— ^iji +r i + 27r 2 r 2 while SUSY technique gives 1 — i?2 as an exponentially decaying function. 

As obvious from eq.(47), R 2 approaches GUE limit as A — > oo but, for finite A- values, it is very different from both 
Poisson as well as GUE. This indicates that the ensembles with distribution parameters giving rise to a finite A do 
not reach stationarity even for infinite size of their matrices, and, their properties being very different from those of 
the equilibrium ensembles, can be referred as "critical". However in our technique, as shown in previous sections, the 
difference between various GG ensembles (within same stationarity limits) manifest itself only in different A-values, 
leaving the functional form of various statistical measures unaffected. Thus RP model as well as ensemble (44), both 
being GGEs and lying on Poisson — > GUE curve, would follow similar formulations for various statistical measures; 
For example, R 2 for both of them is given by eq.(47) although with different formulas for A and both can show the 
critical behavior. However a contradiction arises when one considers the Number variance statistics E 2 (r) which can 
be expressed in terms of R 2 (r) jl2| , 

£ 2 (r; A) = r - 2 f (r - s)(l - R 2 (s)) ds (48) 
Jo 

and therefore should show a similar behavior, as a function of A, for both (RP model and ensemble (44)). But a 
detailed study of RP model by SUSY technique |59| suggests that although it shows critical statistics for fx — cN, 
it can not support linear nature of S 2 = \r with x < 1. As claimed by this study, the difference in S 2 (r) behavior 
arises due to difference in large-r, (A >> r >> vA)j behavior of R 2 (r) in the two cases. 

As our technique is equally well-applicable to both these systems, multifractality should exist in either both or none 
of them. Note that the multifractal nature of an ensemble is so far believed to be indicated by its S 2 -behavior. But 
the latter is not yet clearly understood for RP model (see |3!^,|l|) and therefore question of multifractality is still not 
fully settled. Also note that the earlier results for both models are obtained by SUSY technique using saddle point 
approximation at various stages which may also be misleading, ft is also possible that (i) the seeming multifractality 
of ensemble (44) is the erroneous conclusion of various approximations, or (ii) £ 2 (r) rs x r is n °t always a correct 
indicator of multifractality and therefore its absence in RP model. 

We believe that the £ 2 (r)-behavior is a bigger suspect j3^.4l[]. Our belief has its roots in the direct applicability of 
our technique to Anderson model too. Here also the ensemble for H is located between Poisson — > GUE (for a time- 
reversal breaking disorder) with corresponding i?2-behavior given by eq.(38). Thus for finite A-values corresponding to 
critical region, the eigenvalue statistics is different from Poisson or GUE. But again for E 2 obtained by using eq.(38), 
S 2 (r) x r with x < 1 and therefore if it is indeed an indicator of multifractality of eigenfunctions, our technique 
would suggest its absence in Anderson model. However the existence of multifractality among the eigenfunction of 
Anderson Hamiltonian is experimentally confirmed. 

Our results indicate that multifractality will either be a common feature of all the Gaussian ensembles with finite 
A-values in limit N — ► 00 or it does not exist in any of them. Thus the questions related to critical statistics, the 
correct criteria for multifractality and its analysis by SUSY technique require further probing. 

IV. OTHER CASES 

A. A perturbed Hamiltonian with GG type perturbation 

The intimate connection of RMT — > CM Hamiltonian continues also for system H = Hq + xV with a random per- 
turbation V drawn from a GGE (i.e p(V, y, b) — Cexp(— J2l=i J2k<i a kl;s(Vki; S — bki- s ) 2 )- In this case, the eigenvalues 
evolve due to changing strength of perturbation too. To obtain the desired evolution equation, therefore, one needs 

to consider the sum + £ s T,k<i (7 ~ Vki-s) 2y H;s ^- - b k i- s -^- \ which leads to following equality 



dP 

dx 



s k<l 



dP , dP 



^2^2 (l ~ Vkl; S ) 2y H . s - b kl;s-ZT = ^^hl.s - ^2^2ykl;shl; S (49) 

Oy k l;s O0kl:s] g k<l s k<l 



where Iki- S is still given by same form as eq.(12) but with H replaced by V. As the right hand side of eq.(49) is 
same as that of eq.(14), one again obtains obtains the evolution equation (17) but now Y is given by the condition 



12 



tllclt qy 



OP 

dx 



be given by the following relation: 



Y 



OP 

C>ykl;t 



dP 



Proceeding just as in section I.C, Y can be shown to 



2N 2 + 1) 



EE 



k<l a=l 




c 



(50) 



Again the steady state is achieved for Y 
solution for P is given by Yii<j If 1 



oo which corresponds to x — > oo and yw,s — * 7! the steady state 



Here only x — > oo (with 



dp 



and H = xV) no longer 



represents a steady state, as in the case when V belongs to SGE, but represents a transition state with §y 7^ 0. 
Note from eq.(50) that Y — > 00 as x — > 00, seemingly implying that the equilibrium is reached and therefore H 
belongs to SGE. But, as obvious from H = Hq + xV, in limit x — > 00, if = xV^ and therefore iJ must be a GG 
matrix. This contradiction is a result of the error made in not ensuring the mean spacing of H same as Hq and V 



Here, to ensure the latter, we need to use a modified Hamiltonian, given by H = e « Hq 



-2t 
1-G N 

N 



V with 



-N 1 \n cos(x/N) (same as before in large- iV limit). The effect of this modification on F-P equation (17) is that 



now 



dP 



BP 



dY ~ dr N(1-c- 2t / n ) 

is now replaced by N^ 1 ^ (see eq.(13) of 



dP 



„,-U dP 

7 6«;-3J5^: 



and the coefficient 7 of the drift term 



The Y can now be obtained by the second method given in section I.C. 



B. Non-Gaussian Ensembles 



As mentioned before, the RM models of complex systems can, in general, be non-Gaussian, e.g. p{H) = 

Cexp ^k<i kl w ith / as an arbitrary function and it is not an easy task to obtain the correlations in this case. 
However this case can be analyzed by our method if / is a well behaved function and can be expanded in a Taylor's 
series. To understand this, let us consider an ensemble of real-symmetric matrices H with distribution of a more 
general nature e.g. / as a polynomial of H with degree 2M, fki(x) — Y^l I =ilki{ r )x 2r with C as the normalization 
constant and variances for the diagonal and off-diagonal matrix elements chosen to be arbitrary. 

To obtain an evolution equation in this case, we now consider the sum 2 £ r=1 r J2k<i (7 — Ukl(^)) DM ( r ) dy^TF) 
P = CP and yki(r) = gkijki{r)) where the derivative of P with respect to 7w(l) can be shown to be following (with 
P = Cp) 



dP 



1 



d lk i{l) 27«(1) 



Now as = ~ H llP and mt = - 2 Z A r Ln(r)HZ 



U 6 ^ 



dp _ 



^ TIh(1) 



JY 



£5 dlki{r) 



(51) 



eq.(ll) can be rearranged to show that 2Y^f=i r lki(r) 
(s) on quantities). The required evolution equation in this case, can be obtained from the following equality: 

dP 



dP 

dfkl(r) 



1 p, the second integral in eq.(51) being equal to 

= Iki with Iki given by eq.(12), (but without subscript 



BP 
dlkl (r) ' 



M 



r=l fc<2 



Vki{l)) yu{r) 



dyki (r) 



7E /fei 

k<l 



k<l 



(52) 



where, again, J2k<i hi = £„ wh^P) and 



k<l n 

with Jki now given by following relation: 



dfl n 



d 

dp n 



E 



8 r N 

j -=-e^e/ii^- a *) 



k<l " i=l 



dH k 



M 



k<l 



M 



9ki 



^2(r+l)y k i(r + l] 



dP 
dyki (r) 



M 



(2r + l) + 2^ S y fei (s) 



pdH 

dP 

dyu(s) 



(53) 

(54) 
(55) 
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using these relations as before, one again obtains the F-P equation for P similar to eq.(17) with (3=1 and 
# = 2E fe <*£^iiMr)a!fe where h kl (r) =2r y kl (r) ( 7 - I/w (i)) + ( r + i)(2r + l) y k i(r + l) g H +2(r + l) y fci (r + 

Uki( s ) dy^F) ' Note the condition for V here includes terms of type g y ^f( r ) dy^(s) an< ^ ^ can n0 l° n g er 
be obtained by methods given in section I.C. 



C. Block-Diagonal Ensembles 

The eq.(7) and, therefore, evolution equation (17) of P(/i, Y~) is no longer valid if the matrix H is in a block- 
diagonal form. This is because the eigenvalues belonging to different blocks don't repel each other, are not correlated 
and undergo an evolution independent of the other block. For this case, the evolution of eigenvalues in each block 
can be considered separately, leading to one F-P equation similar to eq.(17) for each block. A detailed discussion of 
this case in given in jl4| . 

V. AN ALTERNATIVE EVOLUTION EQUATION FOR THE EIGENVALUES 

In section LB, the eq.(17) governing the evolution of the eigenvalues was obtained by using the relation (14). 
However, as obvious from eq.(13), P also satisfies the relation 



E 

k<l 



dP 



bkl;. 



dP 



dbki- s 



k<l 



kl:. 



(56) 



and, therefore, one can define a function Z(y k i ;s ,b k i- S ) such that 



(57) 



Here Z is given by the condition ^& 



show that Z = In 



s =l \Vkl;s\b kl . s 



J2k<j 
c 



dP 

dy k i;.. 



?kh. 



dP 



which can be solved (as in section II) to 



rifc<z ris=i 

The eq.(57) also describes the evolution of eigenvalues for the same ensemble (3). But now the "time"-scale is 
such that the eigenvalues seem to be drifting only, hiding the repulsion between them. Again the steady state of 
eq.(57) is given by \Z — Z \ — > oo and the final ensemble as Poisson (with finite, non-zero variances for diagonal 
matrix elements and zero variances for the off-diagonal ones). The ensemble G will now appear as an intermediate 
point in a transition from some initial ensemble — > Poisson ensemble and, in principle, the transition can be used 
for the analysis of G. For example, the critical parameter for Anderson transition (same model as used in section 
III) can be obtained by taking the initial state "O" as metal with energy level distribution described by a GUE 
(< ol >-- 



< Ofj >= 0) and all < O ir , s >= t M which gives Z = ^ [N\na M - 2K\n\t M \] + C The critical 
region will therefore occur as as intermediate point in the GUE — > Poisson transition with transition parameter 



-l 



A = D 2 {Z — Zn) = ttkt A^ln-^ 2Khi T rK . As obvious, the increase of diagonal disorder (-2- < 1) for a fixed 

V u ' 2N [ um 1 1 a / 1 J ' ° y am ' 

hopping rate (t = tj^) will ultimately lead to Poisson statistics, implying localization of states; note here the transition 
occurs backwards in "time" A. However the results for correlations associated with SGE — > Poisson transition are 
not known which leaves eq.(17) as a better tool to analyze the properties of GGEs. . The eq.(17) has one more 
advantage over eq.(58): the reduction of former to CSM Hamiltonian reveals the underlying universality of statistical 
formulation among various complex system. 



VI. CONCLUSION 



In this paper, we have described a new method to analyze the statistical properties of the RM model of complex 
systems. Our technique is based on the exact reduction of spectral analysis in the general case to the one in SGE. 
This greatly reduces the degree of difficulty of the original problem as many of the properties of SGE are already 
known. This also indicates that a thorough knowledge of the properties of SGE or CSM will be highly advantageous 
even for systems with interactions too intricate to be modeled by SGE. So far, the probing of GGE is carried out only 
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by SUSY technique which requires a saddle point approximation at various steps and is not easily applicable, even 
approximately, to cases where our technique can be used for exact probing. Note the main term in GGEs responsible 
for the correspondence with CSM Hamiltonian is due to the repulsion between eigenvalues. As the mathematical 
origin of this term lies in the transformation from matrix space to eigenvalue space which is same for all the hermitian 
ensembles (belonging to same symmetry class) , the correspondence with CSM Hamiltonian should exist for almost all 
of them irrespective of the distribution of their matrix elements. As discussed in section III, our study also confirms 
the conjecture regarding the one parameter scaling of localization and provides the formula for relevant parameter. 

The reduction technique presented here raises some basic questions. Why the reparametrization of the spectral 
properties of different RM ensembles results in to a similar mathematical formulation for them? In other words, why 
the eigenvalues of quantum operators associated with complex systems evolve in a similar ordered way (like equations 
of motion for Calogero particles) notwithstanding the varied nature of their complexity? The reason may lie in the 
following. The eigenvalues and cigcnfunctions of a Hamiltonian evolve due to a change in either degree or nature of its 
complexity. The evolution of the cigcnfunction is chaotic in the sense that the overlapping between the cigcnfunctions, 
associated with two Hamiltonians even with slightly different complexity, decreases rapidly in time (page 2 of Q). 
However an eigenvalue of an operator is its average value in the state described by the associated eigenfunction and 
an ordered evolution of the former will, in general, imply an ordered change in the average behavior of the latter. 
Thus it seems that the eigenvalues and eigenfunctions, on an average, are not able to view the fine subtlities of the 
varied nature of complexity and therefore are not affected too drastically to loose correlations even when nature of 
the complexity changes. Note, for a small change in the interactions, this result is not surprising and used as the base 
for the perturbation theory. But the results in this paper imply that the eigenvalues (and physics based on them) 
even after a violent change in the interactions remain correlated in the parametric space. Thus it seems that certain 
physical properties, based on average behavior of eigenvalues and eigenfunctions, of one complex system are related 
to the physics of another, very different in nature of the interactions. 

I am grateful to B.S.Shastry and N.Kumar for various useful suggestions during the course of this study. I would 
also like to thank K.Frahm and V. Kravatsov for useful criticism of the work. A brief discussion with B. Altschuler 
and B. Huckenstein has also been helpful. 



APPENDIX A: PROOF OF EQS. (3,4,5) 



The use of the eigenvalue equation HU = C/A, with U as a unitary matrix and A the eigenvalue matrix, leads to 
following: 

Y J H ij U jn = X n U in and ^ffy tfo = \ n U* jn (Al) 

3 i 

where i?y = Hij-i + iHij-2- Differentiating both sides of above equation with respect to H k i- S (with s = 1 or 2), we 
get 

^dH kl , s " ^ ,n dH M . a n dH kl , s + dH kl , s U 1 1 

3 3 

Now as '^2 i U* n Ui rn = 5 nm , multiplying both the sides by U* n followed by a summation over all i's, we get the 
following 

t,x " E^T^-^" fA;!i 



1,3 



which further gives 



dHki- s gu 
This can further be used to show that 



" -I s - 1 — [u ln u* kn -{-iyu; n u kn ] (A4) 
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\ " \ " 9X n \ i 

f-^ ^-f differ f-^ fe j 



fc<Z s=l 



k<l 



UlnU* kn J2 i^Hkks + Ut n U kn J2 i'-H-l)' +1 Hkl;, 

s s 

k<l 9kl 

J2 ^-H kl U ln U* kn + ]T ^-Ht k U* kn U ln 
k<i m h>i m 

^ HkiUi n U kn = A„ 



(A5) 

(A6) 

(A7) 
(A8) 



kj 



where eq.(A8) is obtained from eq.(A7) by using Hcrmitian properties of H (i? ; * fe = H k i). By using eq.(A4), One can 
also show that 



E^'E 



d\ n d\ r , 



= EE i2(s_1) ^ - {-i)'UiM [UimU* km (-iyu; m u km ] 



— [u ln u* kn u km u; m + u; n u kn u lm u, 



k<l 



9ki 



kmi 



= Zy^jUlnUknUkmUlm = X UknU km X Ul m U* n - 2S 7 , 



(A9) 

(A10) 
(All) 



where eq.(All) follows from cq.(AlO) by writing J2k<i u in U knUl m Ui m = J2k>i u inU kn Uk m Ui m and the last equality 
in eq.(All) is due to unitary nature of U. 



APPENDIX B: PROOF OF EQ.(6) 

Multiplying both the sides of eq.(A2) by U* m (m ^ n) followed by a summation over all i's, we get the following 



• Jm dH k i- s A„ — A TO 4r ' im dH k 

3 1,3 



lkl;s 



a multiplication of both the sides by U rm followed by a summation over all ra's then gives 

dU rn .H—i 1 U rn 



= ^E 



Qkl , A n — X 7l 



(u* km u ln - (-iyui m u kn ) 



(B2) 



E^E 



d 2 K 
dHl. 



k <l s=l ^kl;a fe < ( 9kl dH kl , 

2 



APPENDIX C: PROOF OF EQ.(7) 



EE^ 1 



dU kr^TT. I 9Uln TJ* I iy+1 9U ln 



U ln + i^Uin + (-l) S+1 iw^U kn + (-1)^1^^; 



dH k i. s dH kky 



dHki. 



k<l s=l 

Now by using eq.(B2) and its complex conjugate in eq.(C2) and by summing over s, we get 

E ski E m^- - 4 E ^ E r4r- + u kn u* kn u lm uu 

k <l s =l OU kl;s fe <, 9kl m A„ A m 

= 4 E E r^ir [^^tWJ 

— — A n — A-- 



A n — A„ 



^2u km u kn 



E^« 



(Cl) 
(C2) 

(C3) 
(C4) 

(C5) 



L fe 
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Now by using the unitary relation J^j Uj m Uj m = 1, one obtains the desired relation (7). 



APPENDIX D: PROOF OF EQ.(15) 



The eq.(12) gives us the following, 

k<l s=l n UfinJ i 

The use of eq.(A8) will further simplify it in following form 



k<l s=l Kl ' s 



pdH 



£ £ = E I II 5 ^ ~ ^npdH 

k<l s=l n ^ n J i 



(Dl) 



(D2) 
(D3) 



APPENDIX E: PROOF OF EQ.(16) 



For each s- value, we have the following relation 

N d — f N 

^ ykl;shl;a = ^ J. 9klUkl-s / \\ S(pi 



K) J^" Hki- s pdH 



k<i 



71=1 r fc<Z 



kl:s 



N 



i=l 



8H-, 



kl:s 



2afej ;s &fej. 



E^E^/n * - ^JL™ + E 

fe<; 



pdH 



^ <9m« 2 7 dH k i-s dHki- s 



(El) 
(E2) 
(E3) 



where eq.(E3) is obtained by using the equality: = —2aki- s {Hki- s — bki- s )p an d JfcZ;s is given by eq.(E9). 

By integrating eq.(E3) further by parts, one obtains 



k<l 8=1 



\ ^ \ ^ T \ \ d \ -\ Qki 

s n k<l 1 



E E £ E f / II *(« -«^ + EE* 

n r k<l 



dH kl 



dH k i- 



a 2 a, 



khs 



E^E^/n*-.) 



EE 

s fc<Z 



fe<( s 
2 ^ ///,/:. ^ ///./: 



E^/n*-*<> 



E- 



k<l s 



V — - V — 



2 £ 



h 

+ E E Jfei ^ 



(E4) 
(E5) 

(E6) 

(E7) 

(E8) 



k<l s 



where Jki- S can be obtained as follows: 
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N d f N OX 
Jm-s = Vki-shi-s ^2 g — / II ^ ~ pdH (E9) 

= -yu-sOu-s / 7TT T pdH (E10) 

J OMki-s 

r N d 

= Vki-sbki.s / TT^ - X i)7TE^ dH ( E11 ) 

J i O tikis 

/ dp dP 
T\S(ni-Xi)— = -yu-sbki-s— (E12) 



where in eq.(A32), the equality = 2aki- s (Hki-s ~ b a . s )p = — d ® p is used. A substitution of eq.(E12) in 

eq.(E8) now leads to the eq.(16). 



APPENDIX F: A GENERAL METHOD TO OBTAIN Y 

Let us consider a transformation of M = 2A^ 2 coordinates {rj} to another set of M coordinates {Yi}, where rj's 
are various coefficients yki-s (total TV 2 ) and bki- s (total iV 2 ). The Yi's should be chosen such that the right hand side 
of the eq.(18), summing over all y^i-sj's an d frfcz ; s's can be rewritten as 

M dP dP dP M dP 

i k<l ' k<l ' j = l J 

where, for our case, gt(rx, .., Tm) — 2 (7 — r^) if is one of the ?/fc2 ;s , and, <7,*(ri, .., Tm) = ~l r i if r i is one of the 
bu-s- 

Now, as we want = ^|r-, with Y\ = Y, this imposes following conditions on the functions Y^s (as can be 

shown by using the theory of partial differentiation) 

dP ,dP8Y t 

l^ = 1^1^9j(ri,r 2: ..,T M ) — — (F2) 



— % v 1 j 



dYi f-fV 'dYi dr 

and therefore 



M dYi 
y29j(,ri,r 2 , ..,r M )-^ 1 = 5u (F3) 



- 1 1 



According to theory of partial differential equations |L5|, the general 

dz 



solution of linear PDE "^f 1 Pi(xx, x 2 , ■ ■, %m)§^: = R is F(ux, 112, .., u n ) = where F is an arbitrary function and 



Ui(x%, x 2 , ■■, x n , Z) = Ci (a constant), i = 1, 2, .., n are independent solutions of the following equation 

dxi dx 2 dxk dxM dZ 



Pi Pi Pk Pm R 



(F4) 



Thus the general solution of eq.(F3) for each Yj is given by a relation Fj(uxj,u 2 j, ..,UMj) — where function Fj 
is arbitrary and Uij{r\ 1 r 2l ..,rM,Yj) = Cij, (i = 1,2, ..,M) (with c^'s as constants) are independent solution of the 
equation 



dri _ dr2 _ chfe dr M _ dlj 

.9i 92 9k 9m Sij 



(F5) 



The above set of equations can be solved for various Yj to obtain Fj. For Y±, we get the relations Y\ — 
^logi j^ 1 = en (i = 1, .., M/2), Y\ + ilog|r 4 | = en (i = 1 + M/2, ,., M), and therefore Fx satisfies the relation 
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F 1 {Y 1 - ^ogj^, ..,Yx - ilog k ^ 7 | ), Yx +7" 1 log|r M / 2 |, ■■, Yx +7 _1 log|r M |) - 0. The function Fx being arbitrary 
here, this relation can also be expressed in the follwing form: 



Yx = 



M 



M/2 , M 



=M/2+l 



C 



(F6) 



where C is another arbitrary function of constants: for example C = C(^log 



7 ^og\r M \, ^logj 



7 MoglrAfl). 



In-7l 



7 MoglrMl, |log 



|r 2 -7| 



2 lw o |r M -i-7l 

Similarly the variables Y"i, i > 1, can be obtained however their knowledge is not required for our analysis. 



APPENDIX G: 



The choice of 7 is based only on the requirement that yu{0) > yki(G) > 7 for all fc, I. Thus 7 can take any value 
such that 7 < xmnyki(G). Let us consider two such possibilities for 7, 7 = 71 and 7 = 72 and try to evaluate properties 
of G on these curves referred as Tl and T2 respectively. Let the value of Y for G on these curves be Yx and Y 2 where 



Y 



Yn 



1 2 
1 = 2~N? EE 



k<l s=l 

! 2 

EE 



1, VM-. 

—In- L 

2 |yfei;s - 



7i I 



dnfyu;, 



2iV 2 



fc<( s=l 

However Yx can also be written as follows 

1 2 
yi = 9^2 E E 



k<l s=X 



2 |yfei; 



1, Vkh 

-In— 

2 \y k i,s~ 



72 1 



71 



— hi&/d :s 
72 



(Gl) 
(G2) 



W, 



72 1 72 



ki.s 



c 



(G3) 



Now as y' kl . s — yu-.s^ 7^ yki-s, this implies that Yi would correspond to a point, different from Y2, on the transition 
curve T 2 and therefore would give properties for the ensemble G different from those given by Y 2 . This conclusion is, 
however, erroneous and is a result of the rescaling applied only to one point Yx on the transition curve Tx- To get 
the right answer, the whole curve Tx should be rescaled which would require a rescaling of the end-points too and 
therefore changed distances on the rescaled curve (call it T[). Thus the point Yx will appear at the same location on 
T{-curve, relative to end-points, where Y2 appears on T2-curve and therefore both will imply the same properties for 
the ensemble G. 
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